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ABSTRACT 

In this work we use gamma-ray burst (GRB) afterglow spectra observed with the VLT/X-shooter spectrograph to measure rest-frame 
extinction in GRB lines-of-sight by modeling the broadband near-infrared (NIR) to X-ray afterglow spectral energy distributions 
(SEDs). Our sample consists of nine Swift GRBs, eight of them belonging to the long-duration and one to the short-duration class. 
Dust is modeled using the average extinction curves of the Milky Way and the two Magellanic Clouds. We derive the rest-frame 
extinction of the entire sample, which fall in the range 0 < Ay ^ 1.2. Moreover, the SMC extinction curve is the preferred extinction 
curve template for the majority of our sample, a result which is in agreement with those commonly observed in GRB lines-of-sights. 
In one analysed case (GRB 120119A), the common extinction curve templates fail to reproduce the observed extinction. To illustrate 
the advantage of using the high-quality X-shooter afterglow SEDs over the photometric SEDs, we repeat the modeling using the 
broadband SEDs with the NIR-to-UV photometric measurements instead of the spectra. The main result is that the spectroscopic 
data, thanks to a combination of excellent resolution and coverage of the blue part of the SED, are more successful in constraining the 
extinction curves and therefore the dust properties in GRB hosts with respect to photometric measurements. In all cases but one the 
extinction curve of one template is preferred over the others. We show that the modeled values of the extinction Ay and the spectral 
slope, obtained through spectroscopic and photometric SED analysis, can differ significantly for individual events, though no apparent 
trend in the differences is observed. Einally we sfress thaf, regardless of the resolution of the optical-to-NIR data, the SED modeling 
gives reliable results only when the fit is performed on a SED covering a broader spectral region (in our case extending to X-rays). 

Key words. Gamma-ray bursts: general - ISM: dust,extinction 


1. Introduction 


Dust plays a central role in the astrophysical processes of inter¬ 
stellar medium and in the formation of stars (e.g., Mathis|1990 


|Draine|[2003| ). Its obscuring effects can introduce large uncer¬ 
tainties to the interpretation of astronomical observations, but at 
the same time offer us the means to study its physical proper¬ 


*Based on observations collected at the European Southern Obser¬ 
vatory, Paranal, Chile, under programs 084.A-0260(B), 085.A-0009(B), 
088.A-0051(B), 089.A-0067(B) and 09LC-0934(C). 


ties, e.g.; the attenuation of light as a function of wavelength, or 
extinction curve, is strongly dependent on the composition and 
size distribution of the dust grains. The origin and properties of 
dust are still poorly known, especially at cosmological distances, 
where a unique probe of dust can be found in gamma-ray bursts 
(GRBs). 


GRBs originate in galaxies at cosmological distances (e.g.. 
Bloom et al.| 1998[ Jakobsson et al. 2012|l. As they are usu¬ 


ally accompanied by bright optical and X-ray afterglow emis¬ 
sion (e.g., |Kann et^L||201l|, they can be used as a powerful 
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tool to study environments at different stages of the Universe’s 


evolution. According to the standard theory (e.g., Gehrels et al. 
2009 1 Gomboc|[2012|l, GRB afterglows are a result of an inter¬ 


action between highly relativistic ejecta, produced in the pro¬ 
genitor’s explosion, and an interstellar medium. The resulting 
emission is of a synchrotron nature, which in its simplest case, 
has a power-law dependence in both time and frequency (e.g.. 
Sari et al.|1998|l. The deviation from a simple power-law is then 


attributed to the absorption and scattering of light by dust grains 
(e.g., |Kann et al.|2006 |201 l| l. GRB afterglow emission is thus 
better suited for studying the extinction, as compared to more 
complex spectra of quasars or galaxies. Modeling of afterglow 
spectral energy distribution (SED) in optical-to-X-ray spectral 
range provides us with information regarding the intrinsic after¬ 
glow spectrum and the properties of the intervening dust. The 
dust properties in random lines-of-sight (LOS) in high-redshift 
galaxies are not known a priori, therefore representative extinc¬ 
tion curves are usually assumed in the modeling. Best studied 
and widely adopted are the average extinction curves observed 
in the Milky Way (MW; 


Cardelli et al. 


Fitz¬ 


patrick & Massa 2007|l, Large Magellanic Cloud (EMC; Pei 


1992| [Gordon et al.||2003|l and Small Magellanic C loud (SMC; 
Prevot et al.||1984[ (Pei||19^ [Gordon et al.||20()^ . The three 


curves differ in their UV slope (steepest in the SMC and shal¬ 
lowest in MW curve) and the prominence of the 2175 A bump 
feature, the latter being the strongest in the MW-type and disap¬ 
pearing in the SMC-type dust. 

Afterglow SED modeling, either limited to the near-infrared 


to ultraviolet (NIR-UV) (e.g.. 

Galama & Wijers 

20011 [Strattal 

et al. [2004[ [Kann et al. 2006 

2010 

Liang & Li 

2010[l or ex- 

tended to the X-ray spectral range ( 

Schady et al. 

2007 2010 


[Greiner et al.|2011) , has revealed systems with mostly low LOS 
extinction where the extinction is preferentially described with 
the SMC-type dust. Studies dedicated to subsamples of more 
extincted aftergl ows ([Kriihler et al.|2011[ [Zafar et al.|2012 [Per-[ 
ley et al.|2013[ Fynbo et al. 2014| l indicated that some SEDs 


show signatures of the 2175 A absorption feature. To get a clear 
and unbiased picture of the extinction properties, [Covino et al.[ 
( 2013[ l analysed the SEDs of a complete sample of GRBs, which 
were not biased towards optically bright events ( [Salvaterra et al.[ 
2012[ Melandri et al.|20l^ . They find that ~ 50% of the after¬ 


glows are found within lines-of sight of low extinction (Ay < 0.4 
mag) and only ~ 13% are heavily extincted (Ay > 2 mag). 

Most of the extinction studies have been done using photo¬ 
metric SEDs. An accurate extinction measurement requires si¬ 
multaneous high-quality data in a broad spectral range from NIR 
to X-rays. Examples of reliable studies of homogeneous data 
sets are those using the observations done with the GR0ND[^ 
instrument ( Greiner et al.|2008 1 , whose capability to simultane¬ 
ously observe in 7 bands in the nIR-UV is especially well suited 
for SED studies ( [Greiner et al.|20lT[ [Kriihler et al.|201 l| l. 

Extinction studies could be improved by using spectroscopic 
instead of photometric SEDs. Zafar et ^ ([201 1[) (hereafter 


Zll) studied a sarrmle of 41 optical afterglow spectra obtained 
mostly with FORl^ (VET). They showed that the SEDs with 
prominent 2175 A absorption are among the most extincted with 
Ay > 1.0, while the events with SMC-type dust have low ex¬ 
tinction (Ay < 0.65). This led them to conclude that the low 
detection rate of afterglows with EMC or MW-type dust is more 
likely due to observational bias against dusty LOS than due to the 
MW-type dust to be rare in high redshift environments. While 


of a limited spectral coverage, the spectroscopic SEDs allowed 
them to prove convincingly that the afterglow spectra are indeed 
consistent with a simple synchrotron model. 

In this work we take the next step in the SED modeling by us¬ 
ing the afterglow spectra obtained with the X-shooter instrument 
( [Vernet et al.[|20ir| ). X-shooter is a state-of-the-art intermedi¬ 
ate resolution spectrograph, mounted on the VET, that simulta¬ 
neously covers a broad spectral range with three spectroscopic 
arms: ultraviolet (UVB; -3000-5500 A), visible (VIS; -5500- 
10000 A) and near-infrared (NIR; -10000-25000 A). The large 
spectral coverage of X-shooter spectra offers us a unique oppor¬ 
tunity to apply a detailed extinction curve analysis. The power 
of X-shooter in extinction curve studies was illustrated by the 
recent observation of GRB 140506A, whose spectrum revealed 
a unique extinction signature in an afterglow spectrum: a very 
strong flux drop below - 4000 A (in the GRB’s rest system) is 
unprecedented in the study of GRB environment and has been 
found only in a few other LOS to other types of objects so far 
( [Fynbo et al.[2014 1 . Our aim is to use afterglow spectra acquired 
with the X-shooter instrument to derive the dust properties in 
the LOSs of GRB host galaxies and to evaluate the applicability 
of the commonly used extinction curves in this type of analysis. 
The sample, data preparation and analysis are presented in Sec¬ 
tion!^ Results and detailed discussion are given in Sections 
and W respectively. We summarize our conclusions in Section!^ 


Throughout the paper the convention Fy(t) 


-a^-13 


adopted, where a and p are temporal and spectral afterglow 
slopes, respectively. Times are given with respect to GRB trigger 
time. 


2. Data and Analysis 

2.1. Preparation of X-shooter and X-ray spectra 

As a part of the X-shooter GRB GTO prograrr0 spectra for - 60 
GRB afterglows have been acquired in the period between 2009 
and 2014. From this sample we selected those GRBs whose 
spectra, according to photometric observations of the afterglows, 
are not contaminated by supernova or host galaxy emission. The 
spectra were reduced and calibrated using version 2.0 of the X- 
shooter data reduction pipeline ( Goldoni et al.|200'^ Modigliani 
[et al.|2010l ) - details are described in Fynbo et al. (in prep). In 
particular, the instrument’s response function, required to flux- 
calibrate a spectrum, was obtained by comparing an observed 
spectrum of a spectrophotometric standard star with the tabu¬ 
lated values. Flux calibration of the X-shooter spectra has to be 
robust for our science case. The observations of the standard 
stars used for the flux calibration are performed only once per 
night with a 5" wide slit and binning different than the one used 
for the science spectra. Afterglow observations are done with 
much narrower slits, usually with 1.0", 0.9" and 0.9" for UVB, 
VIS and NIR spectra, respectively. Even if the sky conditions 
are the same when the standard star and afterglow are observed, 
a narrow slit loses much more light in the case where seeing is 
comparable to or larger than the slit width, which can be espe¬ 
cially problematic due to the weak wavelength dependence of the 
seeing ( [Roddier|1981[ ). Other potential problems are the loss of 
light in the case the slit is not positioned in the direction of the 
parallactic angle, and a faulty performance of the atmospheric 
dispersion correctors in the UVB and VIS spectroscopic arms, 
which have been disabled in Aug 2012. 


'http://www.mpe.mpg.de/~jcg/GROND/ 

^http://www.eso.org/sci/facilities/paranal/instruments/fors.html 


^PIs: J. Fynbo and L. Kaper 
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Table 1. Presentation of the sample 


GRB 

Z 


N° 

^*H,X 

[lO^^cm-^] 

Tmid 

days 

Arx 

ks 

100219A 

4.667 

0.208 

6.7 

0.55 

15-50 

100418A 

0.624 

0.200 

4.8 

1.47 

80 - 250 

100814A 

1.44 

0.054 

1.6 

4.1 

200 - 400 

100901A 

1.408 

0.270 

13 

2.75 

160 - 500 

120119A 

1.728 

0.295 

1.1 

0.074 

4-20 

120815A 

2.358 

0.320 

8.4 

0.086 

6-20 

130427A 

0.34 

0.055 

1.8 

0.70 

40-70 

130603B'*> 

0.357 

0.063 

2.0 

0.35 

15-60 

130606A 

5.913 

0.066 

2.0 

0.329 

10-40 


Notes. GRBs in the sample. For each GRB we report the basic informa¬ 
tion required in the analysis: redshift, Galactic extinction A® and equiv¬ 
alent neutral hydrogen column density in the burst’s line-of-sight, 
mid-time of the X-shooter observation Tmid and the time interval ATx 
(in ks) used for the construction of the X-ray part of the SED. Times 
are given relative to the start of the GRB y-ray emission in the observer 
frame. GRB 130603B belongs to the short class, while other GRBs in 
the sample are of the long class. 


For these reasons we cannot blindly rely on the instrument’s 
response function obtained in the calibration process and we 
therefore require multiwavelength photometric data of the stud¬ 
ied objects to check and validate the flux calibration (see also 
Kriihler et al. in prep. ). Our sample is thus limited by the avail¬ 
ability of multi-band photometric data in the literature. Photo¬ 
metric observations should be available at or near the epoch of 
X-shooter spectra. We found that for several GRBs in the full X- 
shooter sample not only the absolute flux level, but also the flux 
calibration as a function of wavelength did not match the pho¬ 
tometric SED. Since it is extremely difficult to reliably account 
and correct for all the effects (either of observational or techni¬ 
cal nature) which influence the calibration, we decided to work 
only with the spectra with flux calibration where, after apply¬ 
ing a correction to the absolute flux calibration, the difference 
between the spectrum and photometry is less than 10% in all 
bands. The photometric data are used only to check the validity 
of the flux calibration - we do not use them to correct the slopes 
of the flux-calibrated spectra. UV photometric observations are 
seldom available around the X-shooter epoch. The seeing for all 
but one event, as measured from the 2D trace of the spectra, is 
small enough that the slit losses as a function of wavelength at 
the blue part of the UVB arm should be negligible and there¬ 
fore we trust the calibration of the bluest SED part. The seeing 
was bad during the observation of the GRB 100901A afterglow 
{FWHM ~ 1.5 ), but in this case the UV photometric observa¬ 
tions were available, confirming the calibration in the UV part of 
the SED is fine. In summary, a GRB was included in our sample 
if: 

1. multwavelength photometric observations of its afterglow 
around the epoch of X-shooter observation are available, 

2. its spectrum is not contaminated by host or supernova emis¬ 
sion, 

3. the difference between the flux-calibrated spectrum and pho¬ 
tometry is less than 10% in all bands. 


Our final sample consists of nine GRBs, eight of them belonging 
to the long and one to the short clas^ The sample is presented 
in Table [T] In the future, more published light curve data will 
enable us to expand the analysis to a bigger sample. 


The absolute flux calibration of the spectra in this sample 
was first fine-tuned with the photometric observations. Spectra 
were corrected for the extinction originating in our Galaxy us¬ 
ing ICardelli et al.| ( |1989| l extinction curve (assuming the ratio 
of total-to-selective extinction Ry = 3.1) and Galactic extinc¬ 
tion maps ( Schlafly & Einkbeiner|20lT| ). Regions of telluric ab¬ 
sorption and strong absorption lines originating in the GRB host 
galaxies were masked out. Spectra were then rebinned in bins 


of widths 30 - 100 A in order to reduce the noise and to guar¬ 
antee a comparable weight of the optical and X-ray SED part. 
The binning was performed with the sigma-clipping algorithm; 
all data points within a binned interval that differ for more than 
three standard deviations from the mean of the points in the in¬ 
terval were rejected. Since the flux in the binning region is not 
constant (i.e., the spectra follow power-laws), we paid special 
attention to avoid removing the tails of the binned region during 
the sigma-clipping. Due to a large number of points that were 
binned we calculated the errors as an average of a standard devi¬ 


ation of binned points from the mean value cr/ ViV, where N is 
the number of binned data points. We cross-checked these val¬ 
ues with a Monte Carlo simulation. We assumed the errors of 
the data points, as obtained in the reduction and calibration pro¬ 
cedure, are Gaussian. We then resampled the data points in each 
interval for a thousand times and computed the l-cr equivalent of 
the resulting distribution of mean values. The errors, computed 
by the two methods, are comparable. 


The reddest part (e.g., K band) of the NIR spectrum is sel¬ 
dom accurately flux calibrated. This is due mainly to the strong 
vignetting of the K band (section 2.4.9 of the X-shooter User 
Manual) which prevents a reliable sky subtraction especially for 
long exposure and faint sources. In addition, many observations 
were conducted with a special X-blocking fllte]|^ Eor these rea¬ 
sons, only A < 20000 A part of the NIR spectra have been used 
in the analysis. 


X-ray data from the Swift/XRT instrument ( Burrows et al.| 
|2005| l were taken from the online repository of X-ray afterglow 
spectra ( |Evans et al.|2009l l. The X-ray SED was built from the 
light curve in a time interval around the epoch at which the X- 
shooter spectrum was taken: only a portion of the X-ray light 
curve without significant spectral evolution was considered. The 
mean time of the X-ray SED was computed as 2;(f/^h)/ 
where f, is the mid-time of individual exposure and Af, is the ex¬ 
posure length. The light curve at the considered time interval 
was fitted with a power-law function: by knowing the slope, the 
X-ray SED was normalized by interpolation to the epoch of X- 
shooter observations. The uncertainty of calculated normaliza¬ 
tion is never greater ~ 10%; by varying the normalization within 
this uncertainty for each analysed SED we found that the best 
fit parameters are consistent within errors. X-ray channels were 
rebinned in order to have at least 20 photons in each channel. 


'^Traditionally, GRBs are classified into a long or short class accord¬ 
ing to their observed duration (i.e., longer or shorter than ~ 2 s) and 
spectral properties jKouveliotou et al.|1993| l. However, a reliable clas¬ 
sification into the two classes, which correspond to different progenitor 
types, is usually more complicated (e.g., Zhang et al. 120091. 

^The filter blocks the spectral range above 2 pm and thus prevents 
the scattered light from X-band orders to contaminate the J- and H- 
band background (e.g., Vernet et al. 2011 1 . 
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which ensures that the data are roughly Gaussian and allows a 
reliable use of statistics. 


2.2. SED modeling 

In most cases the optical-to-X-ray SED of the afterglow is stud¬ 
ied hours to days after the GRB explosion (see Table [^l, when 
the shape is expected to depend on the relative values of the 
cooling (Vc), optical (vo), and X-ray (vx) frequencies ( |Sari et ah 
1998| l. The SED can therefore be described with a power-law 
function Fy - (i.e., Vc < vo or Vc > vx) or, in case 

Vo < Vc < Vx, by a broken power-law; 


Fy = E(, 


v-Po 


V < Vc 


-/3y s,. 

Vc V y > Vc, 


( 1 ) 


where Afi - fix- Po- 

The optical afterglow is attenuated due to light scattering by 
dust particles and can be accounted for in our observations by: 
Ev.obs = FylO^^-^^', where is a wavelength-dependent ex¬ 
tinction in the host galaxy frame. We assume the extinction is 
dominated by the dust in GRB hosts rather than being attenu¬ 
ated by intervening galaxies that happen to occur in the LOS be¬ 
tween the GRB event and observers on Earth - there are known 
cases when other galaxies are found in the LOS (e.g., |Sparre| 
et al.|2014| l, however in most cases the host galaxy is the dom¬ 


inant absorption system (e.g., |Vergani et al.|2009[|Schady et ah] 
|2010| l. The value of - Ayf{A',Ry, ...) is parametrised us¬ 
ing the rest-frame extinction Ay (extinction in the rest frame V- 
band, at 5500 A) and the extinction curve /, which reflects the 
properties of dust in the line-of-sight. We consider the average 
observed extinction curves of Milky Way (MW) and Large and 
Small Magellanic Clouds (EMC, SMC; |Pei|1992) in the analy¬ 
sis. Initially we also considered the average extinction curve of 
starburst galaxies ( |Calzetti et al.|2000| ), but we find the starburst 
curve is inadequate in describing extinction in lines-of-sight of 
GRB hosts - in addition to being statistically unjustified, fits of 
starburst models for all events in our sample result in physically 
unreasonable best-fit parameters of p and Ay. The conclusion 
that starburst models cannot describe dust in GRB lines-of-sight 
was also reached in other studies (e.g., |Covino et al.|2013 1 . This 
is not surprising as it is not expected that extinction in LOS to in¬ 
dividual sources and galaxy integrated extinction properties are 
the same. 

Lining only the X-shooter SED may result in a degenerate 
solution between the values of the spectral slope, p, and the vi¬ 
sual extinction. Ay, i.e., data could be successfully represented 
both by a shallow slope and high amount of extinction or a steep 
slope and low extinction. The problem can be resolved by con¬ 
straining one of the parameters. We decided to consttain the 
spectral slope by including the X-ray part of the SED in the 
fit. By analysing a large afterglow sample, Z11 found that the 
difference between the optical and X-ray slope is very close to 
the theoretically predicted value of AjS = 0.5 in the majority of 
cases. Nevertheless, we left the slope difference A/? as a free 
parameter of our model. The photoelectric absorption of soft X- 
rays by metals is assumed to originate in the Galaxy and in the 
host galaxy at a known redshift z. Galactic equivalent neutral 
hydrogen column density taken from maps provided by 


Kalberla et al. 


( 2005[ p| while the value in the host galaxy line-of- 


sight Ah.x is left as a free parameter of the model. Eollowing the 


®In general, in addition to the atomic hydrogen Hi, provided hy 
maps of |Kalberla et al.|(2005^ , a contribution of molecular hydrogen H 2 
should be taken into account when estimating Galactic X-ray absorption 
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discussion of 


Watson et al. (|2013| l we assume solar abundances 


o: 


f |Anders & Grevesse| ( |1989| l. 

In principle, we do not know whether the spectral break is 
sharp as in Equation [T| or t he transition is more mild and smooth 
(e.g., Granot & Sari 2002[ Uhm & ZhangpOf^ . Indeed, at least 


in one case the cooling break appears to be smooth: Eilgas et al. 
P011| l found that multiepoch SEDs of GRB 091127 can be de¬ 
scribed with a broken power-law with a smooth spectral tran¬ 
sition. We carefully evaluated the effect of the smoothness on 
our modeling of broadband SEDs with data of X-shooter qual¬ 
ity. We find that the available data cannot constrain the magni¬ 
tude of the smoothness (i.e., the sharpness index in Equation 1 
in |Eilgas et al.|20lT] l. Eurthermore, if a break is very smooth and 
lies near the optical region, the values of other parameters (like 
Ay and Pq) obtained in the modeling are considerably affected 
and uncertain. This is not the case if the break is sharp. In the 
absence of a better knowledge on the spectral smoothness, we 
assume that the breaks are sharp in order to reduce the possible 
systematic parameter eiTors. 

Eor GRBs occurring at z > 2, the host’s Lycr absorption line 
enters the X-shooter observational window. In order to better 
constrain the UV slope of the extinction curve we decided to in¬ 
clude the modeling of the red wing of the Lya line in our analy¬ 
sis. The Lya line is characterized by its central wavelength Ai^ya, 
column density of the absorbing gas Ahi and Doppler parameter 
b. As expected for the damped Lyman alpha absorbers, the fitting 
turned out to be very insensitive to the latter, therefore we fixed 
its value to b - 12.6 km s \ which is the average Doppler pa¬ 
rameter of GRB host galaxy absorption lines (|de Ugarte Postigo 


|et al.|2012) l. To model the Lya line we used the analytical ap¬ 
proximation derived by |Tepper-Garcfa| ( |2006| |7007| l. The final 
model applied to the data and already corrected for Galactic ex¬ 
tinction and photoelectric absorption can therefore be summa¬ 
rized as: 


Ev.obs = Ey X 10 X exp[-AH,xO'(v)] x exp[-TLya(v)], (2) 

where cr(v) is a cross-section for photoelectric absorption oc- 
cmi'ing from the gas in the host galaxy. 

The SED fitting was carried out with the spectral fitting pack¬ 
age XSPECvl2.8 ( |Arnaud|[1996) l.The analysis was done sepa¬ 
rately for (/) the X-shooter spectrum and (ii) the broadband SED 
including the X-shooter and X-ray spectra. We also modelled a 
broadband SED with photometric data instead of the X-shooter 
spectrum, in order to compare the two types of analyses. Confi¬ 
dence intervals were computed at 90% confidence level follow¬ 
ing |Avni| ( |1976[ ) and |Cash| ( |1976| l with one parameter of inter¬ 
est. Confidence intervals were being computed independently 
for each model parameter. We did not investigate possible corre¬ 
lations of uncertainties between the model parameters. We con¬ 
sidered a fit to be successful if the null probability (that is, the 
probability of getting a value of as large or larger than ob¬ 
served if the model is correct) was better than 10%. If the broken 
power-law provided a better fit than a single power-law, we used 
the E-test (with the probability of 5% as a threshold) to assess 
whether the improvement in the;^^^ is statistically significant. 


( [Willingale et al.|2013^ . We ignore the molecular component in order 
to make our modeling of the X-ray spectra comparable to analyses of 
larger samples in previous studies (e.g., see Section [4^ and Figure]^, 
in which the H 2 has not been considered. Nevertheless, we checked that 
taking the simplified model does not affect our global modeling, i.e., the 
values of spectral slopes and extinction. 
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!^rest[10'^Hz] [10l= Hz] 


Fig. 1. Rest-frame afterglow SED analysis of GRB 120815A. Broadband SED is shown in the left plot. Optical and X-ray data are plotted with 
black points, the best fitted model is plotted with solid blue lines and the intrinsic afterglow spectrum with dashed blue line. Zoomed part of the 
X-shooter SED is shown in the right plot, where residuals to the best-fitted model are also shown. Vertical lines divide the spectral regions covered 
by the NIR, VIS and UVB spectrograph arm. Orange points are photometric observations, used to calibrate the absolute flux of the spectrum. 
References for photometric data points and details on the fitting procedure are reported in Section]^ Similar plots for other GRBs in our sample 
are shown in Figures |^and|^ 


3. Results 

Best-fitted model parameters for the SEDs of the GRBs in our 
sample are reported in Table An example of the best model 
in the case of GRB 120815A is plotted in Figure plots corre¬ 
sponding to other GRBs are shown in Figures]^ andDetailed 
fitting results are outlined in Table[AT]and discussed in this Sec¬ 
tion on a case by case basis. 

In the following we report qualitative fitting information for 
each analysed SED in the sample. We also briefly discuss the 
compatibility of our SED results with the closure relations for 
each case, assuming models assembled by |Racusin et al.| ( |2009[ 
Table 1). To simplify the discussion we do not consider struc¬ 
tured jet models and models with continuous energy injection. 
We caution that for a detailed understanding of each event, the 
reader should consult works dedicated to each GRB. 


3.1. GRB 100219A 


Photometric data used for normalization and photometric SED 
analysis were obtained from Thone et al. ( |2013 1. The X-shooter 
SED alone can be fitted by a single power-law and SMC or EMC 
extinction curve. The fit is improved when we include X-ray 
data - the broadband SED can be fitted by a power-law and all 
three extinction curves. The best fit is achieved with the EMC 
extinction curve. A broken power-law does not statistically im¬ 
prove the fit, therefore we do not consider it further. This is 
a high-redshift GRB (z = 4.667), therefore we had to include 
the Lya absorption to constrain the UV slope. The value of the 
log Ahi = 21.0+0.1 that we obtain in the fit is comparable within 
errors with the one derived by [Thone et al. ( |2013| ) from normal¬ 
ized spectrum (i.e., log Ahi = 21.14 + 0.15). Fixing the value to 
log Ahi = 21.14 results in a bit worsestatistics, but otherwise 
the other parameter values do not change. 

For the photometric SED we use only filters not affected by 
host’s Lya and Lyman forest. According to the statistics, broad¬ 
band SED with photometric optical points is best-fitted with a 
power-law and MW curve. However, visual inspection shows 
that the optical SED part is better described with a power-law 
and SMC or EMC extinction. The latter two give a comparable 


power-law slope and slightly lower extinction Ay as the one ob¬ 
tained with the X-shooter broadband SED. A rather low reduced 
values are due to small number of data points (both at optical 
and X-rays) included in the fit. 


Optical light curve is characterized by several moderate re- 
brightenings ( |Mao et ah 2012| l. The X-shooter spectrum is taken 
right after a bump peaking at ~ 20 ks. The light curve slope at 
the X-shooter epoch is a ^ 1.31 ( Thone et al.|2013| . The case of 
an ISM environment and Vc > vx seems to describe the optical 
light curve the best (a ~ 1.1, p ~ 2.5). At this point the X-ray 
light curve is in a transition from a shallow to steep (a ~ 2.9) 
phase, which is not seen in the optical ( [Thone et al.||2013] l. If 
this is due to a geometrical effect (i.e., a jet break) rather than 
spectral evolution, our result is wrong. Due to a combination of 
aforementioned bump and sparse data it is hard to estimate the 
correction to normalization of the X-ray spectrum due to a jet 
break. We approximately estimate that X-ray afterglow at the X- 
shooter epoch would be a factor of 2 brighter in the absence of a 
jet break. After applying this correction and repeating the fit, we 
find that a single power-law and the LMC-type extinction curve 
still describe the data the best, though the extinction (Ay ^ 0.26) 
and spectral slope (p ^ 0.66) do change a bit, as expected. 


3.2. GRB 100418A 

Photometric data used for normalization and photometric SED 
analysis were obtained from de Ugarte Postigo et al. (in prep). 
Three afterglow spectra were obtained at ~ 0.4,1.4 and 2.4 days 
after the burst. Due to the bad flux calibration of the first epoch 
spectra and possible contamination from host galaxy and super¬ 
nova emission of the third epoch of observation we only use the 
second epoch spectrum. The X-shooter SED alone can be fitted 
by a single power-law. However, the broadband SED requires 
a broken power-law shape. The difference in spectral slopes is 
A/3 - 0.31 ± 0.09. Fixing the spectral difference to A;6 = 0.5 
significantly changes only /3o, while other parameter values stay 
almost unchanged - the E-test probability (< 1%) confirms that 
the model in which both spectral slopes are left free to vary is 
statistically better. While the x^ statistics is similar for all three 
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Table 2. Results of the best-fit models of the broadband SEDs 


GRB 

Extinction 

A^h,x 

Pi 

Pi 

{a) 

^break 

log Am (Ly a) 

Av 

A^dof 


type 

[lO^^cm-Z] 



[10'5 Hz] 

log cm“^ 



100219A 

EMC 

<6.2 

0 73+0-02 



21.0!“;J 

0 23+^-^^ 
^•^*^-0.02 

34.8/30 

100418A 

SMC 

-0.12 

0 73+0-07 

1 04+^-^^ 
^•^^-0.03 

0 6+^-^ 
^■^-0.1 


-0.02 

20.8/23 

100814A 

SMC 


0 52+0-07 
-0.07 

1 05 

^•^'^-0.02 

9 9+3.8 
^■■^-0.1 


^•^^-0.03 

70.8/66 

100901A 

SMC 

^•^'^-0.18 

-0.04 

1-06!°:“ 

c q+8.8 
^■^-3.2 


^•^^-0.03 

44.2/41 

120119A 

EMC 

1 98+0.50 

0.89!“;“! 



93 4+0.2 

1 07+003 
-0.03 

106.0/81 

120815A 

SMC 

0.66!°;^ 

0-38!“:“ 

0.84!“;“ 

1 4+0.7 
-0.8 

99 3+0.2 
^^•■^-0.2 

0 32+0-02 
^•*^^-0.02 

26.0/47 

130427A 

SMC 

0.08!“;“ 

0 37+0-05 
-0.04 

0.68!“;“! 

o-7!Si 


0 16+0-02 
^•^^-0.02 

129.3/147 

130603B 

SMC 


Q 42+0-12 
-0.22 

0 Q9+0.08 

^•^^-0.04 



1 1 q+0.23 
^•^^-0.12 

21.3/23 

130606A(*) 

/ 

<3.5 

0 96+^-®^ 
^■^^-0.02 



1 Q Q+0.2 
^^•^-0.2 

<0.01 

48.8/31 


Notes. Detailed summary of the fitting results is outlined in Table lA.lj and discussed in Section|^ 

(a) Host rest-frame value. 

(b) GRB 130606A is found to be consistent with Ay ~ 0, therefore no extinction curve is needed for modeling. 


extinction types, visual inspection reveals that the SMC extinc¬ 
tion curve provides the best fit to the blue part of the optical data. 

The broadband SED with photometric optical points is fitted 
equally well with all three extinction curves and a broken power- 
law. The best-fit slopes are comparable to the values obtained 
from the X-shooter broadband SED fit, while extinction values 
are slightly higher (but within errors). The broadband fit with a 
single power-law also describes the data well, however the E-test 
probability (< 5%) confirms the broken power-law improves the 
fit significantly. 

Initially shallow evolution of the optical light curve was fol¬ 
lowed by a rebrightening, reaching its peak brightness at ~ 0.6 
days ( [Marshall et al.|20lT| ). The data from de Ugarte Postigo et 
al. (in prep) suggest an optical late-time steepening of ao ~ 1.5, 
similar to the X-ray decay of ax ~ 1-4 in this late phase. 
[Marshall et al.j ( [201 l| l find a shallower optical steepening of 
ao ~ 1.1. In the latter case the difference between the opti¬ 
cal and X-ray light curve slopes cannot be explained within the 
standard model without a spectral break between the optical and 
X-ray region. On the other hand, the case of ao ~ 1 5 could be 
explained within a model of non-spreading uniform jet in a wind 
environment. We do not find any significant color evolution in 
a few days around the X-shooter epoch (de Ugarte Postigo et al. 
in prep), which is in contradiction with the stellar wind environ¬ 
ment scenario. 


3.3. GRB 100814A 


Photometric data used for normalization and photometric SED 
analysis were obtained from |Nardini et ak (2014l. Three after¬ 
glow spectra were obtained at ~ 0.038,0.089 and 4.1 days after 
the burst. Due to the discrepancy between flux calibrated spec¬ 
tra and photometric measurements at the first two epochs, only 
the spectrum taken 4.1 days after the burst has been analysed. 
The X-shooter SED alone can be fitted with a SMC extinction 
curve and a power-law intrinsic behavior. Statistically there is 
no need for a spectral break. EMC and MW curves fail to repro¬ 
duce the data. The broadband SED is fitted well with a broken 
power-law and the SMC or EMC extinction curve. However, the 
EMC clearly overpredicts the 2175 A bump. We therefore prefer 


Article number, page 6 of 17 


the SMC extinction curve as the best model to describe the data. 
The difference in spectral slopes is A/3 - 0.53 + 0.08. Eixing 
the spectral difference to A/3 = 0.5 does not change the results 
significantly. 

The broadband SED with photometric optical points is fitted 
equally well with all three extinction curves and a broken power- 
law. The position of the spectral break is similar to the position 
obtained with the X-shooter broadband SED, while the spectral 
slopes are different. The broadband fit with a single power-law 
is bad (yVdof > 3.0). 

Optical afterglow light curve of GRB 100814A is character¬ 
ized by a strong rebrightening at ~ 20 ks. At the X-shooter epoch 
the optical and X-ray light curves decay as a =» 2.25,2.30, re¬ 
spectively ( [Nardini et al.[2014j l. Our SED results are marginally 
consistent with the model of a spreading uniform jet. Light curve 
evolution in optical prior to the X-shooter epoch is chromatic - 
afterglow is becoming redder with time ( [Nardini et al.|2014) l. If 
the position of the cooling frequency is indeed between the op¬ 
tical and X-ray (as we find it to be), then this would suggest that 
the frequency is traveling towards higher frequencies. 

3.4. GRB 100901A 

Photometric data used for normalization and photometric SED 
analysis were obtained from Gomboc et al. (in prep). The X- 
shooter SED alone can be fitted by a single power-law spectrum. 
Broadband fit reveals the need for a spectral break between the 
optical and X-ray regions. SMC extinction curve provides an ex¬ 
cellent fit to the data, while the EMC and MW curves are com¬ 
pletely inadequate to describe the SED. The difference in spec¬ 
tral slopes is A/3 - 0.56 + 0.07. Eixing the spectral difference to 
A/3 - 0.5 does not change the results significantly. 

The broadband SED with photometric optical points is fitted 
best with a broken power-law. Model with the SMC-type dust 
provides a marginally better fit than EMC- or MW-type dust. 
The values differ quite a lot from the ones obtained from X- 
shooter broadband SED, owing to the unavailable near infrared 
photometric observations. 

The afterglow exhibits an extreme rebrightening phase at op¬ 
tical wavelengths (e.g., [Hartoog et al.[[20T3l Gomboc et al. in 
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prep.). Light curve at the time of the X-shooter observation 
evolves practically achromatically with ao = 1.52 + 0.05 and 
ax = 1.55 + 0.05. None of the closure relations can simultane¬ 
ously describe the observed spectral and temporal properties at 
late times. 


3.5. GRB120119A 


Photometric data for normalization and photometric SED anal¬ 
ysis were obtained from [Morgan et aL] ( 2014| l. Two afterglow 
spectra were obtained at ~ 0.074 and 0.20 days after the burst. 
Due to poor signal-to-noise of the second epoch spectrum we 
only use the one obtained at 0.074 days. This is the only case 
in our sample with a clear 2175 A bump in the SED. In addi¬ 
tion, the blue part of the SED is already affected by the red wing 
of the Lya line. We cannot obtain a very good fit modeling the 
X-shooter SED alone. The broadband SED is fitted best with a 
power-law and EMC extinction curve. Even so, the EMC curve 
overpredicts the strength of the 2175 A bump. A broken power- 
law on a broadband dataset (and X-shooter data alone) results in 
an unphysical result of yS 2 < Lya line lies in the bluest part 
of the spectra, where signal-to-noise is very low. This is prob¬ 
ably the main factor contributing to the overestimated value of 
log Ahi = 23.4 + 0.2 with respect to the value measured from the 
normalized spectra (e.g., log Ahi = 22.5 + 0.3; Vreeswijk et al. 
in prep.). 


The broadband SED with photometric optical points is fitted 
best with a single power-law and EMC extinction curve. Best-fit 
spectral slope is similar to the one obtained with the X-shooter 
broadband fit, while the extinction Ay is slightly lower. Statis¬ 
tically there is no need for a spectral break in the fitted spectral 
region. 


The early time light curve is characterized by color evolu¬ 
tion, possibly explained as a result of dust destruction ([Morgan 
[et al.|2014[ ). At the X-shooter epoch the light curve evolution is 
achromatic in the optical-to-X-ray spectral range with a ~ 1.3. 
Spectral and temporal indices are consistent with the case of 
vq > yx and an ISM environment. 


3.6. GRB 120815A 


Photometric data used for normalization and photometric SED 
analysis were obtained from [Kruhler et al. ( 2013| l. The X- 
shooter SED can be fitted with a broken power-law and SMC 
extinction curve. Models with EMC and MW extinction curves 
are successful in reproducing the data, but the values of the post¬ 
break slopes are unrealistically steep (given the X-ray part of 
the SED). The broadband SED can only be fitted with a broken 
power-law and an SMC extinction curve. This GRB originates 
at z = 2.358, therefore we had to include the Lya absorption to 
constrain the UV slope. The value of the logAni = 22.3 + 0.2 
that we obtain is slightly larger than the one derived by Kruhler 
[et al.|P013[l, who found logAni = 21.95 + 0.10. As discussed 


by Kriihler et al.| ( [2013[ l, the Lya line is contaminated by vi- 
brationally excited lines that form a continuum around this 
spectral region and cause the line to appear stronger. Not taking 
the Hj into account, the column density value of (contaminated) 
Lya is log Ahi ~ 22.1, closer to our value. Eixing the value to 
log Ahi = 22.1 in the modeling does not change the results (see 
Table I A. ll . The difference in spectral slopes is Aj3 - 0.46 + 0.07; 
fixing it to A/3 = 0.5 does not change the results of the fitting sig¬ 
nificantly. 


We did not use the g'-band photometric point in the photo¬ 
metric SED analysis because the measurement is affected by Lya 
absorption. Broadband SED with photometric optical points is 
fitted best with a broken power-law and SMC or LMC curve. Eit 
with a MW curve, on the other hand, is successful with a single 
power-law spectrum. In the latter case the visual inspection re¬ 
veals the model does not describe optical part of the SED very 
well. Best fit host extinction is lower than the one obtained from 
the X-shooter broadband fit. 

The light curve of GRB 120815A exhibits a smooth tran¬ 
sition from shallow (a = 0.52 + 0.01) to somewhat steeper 
(a = 0.86 ± 0.03) decay at a; 0.05 days after the burst ( Krtih-[ 
[ler et al.|2013| l. Its behavior is achromatic in the NIR-to-X-rays, 
implying an absence of spectral break in this spectral regime. 


3.7. GRB 130427A 


Photometric data used for normalization and photometric SED 
analysis were obtained from [Perley et al.[ ( [2014[ ). The X-shooter 
SED alone can be fitted by a broken power-law and all three 
extinction curves. The broadband SED also requires a broken 
power-law shape. While the reduced implies the SMC or 
LMC curve provide a fit of a similar quality, a visual inspection 
shows that the model with the LMC-type dust fails to reproduce 
the data in the bluest X-shooter region. We thus prefer the SMC 
curve with low host extinction as the case best describing the 
real conditions. The break is required to occur within the ob¬ 
servable X-shooter range. The difference in spectral slopes is 
Ay6 = 0.31 ± 0.05. Fixing the spectral difference to Af3 - 0.5 re¬ 
sults in statistically worse, yet still acceptable fit. Visual inspec¬ 
tion reveals that the X-shooter continuum in the latter case is not 
fitted that well. This fit also results in a shallower optical slope 
(jSo ~ 0.2) and consequently higher extinction (Ay ~ 0.35). 

The broadband SED with photometric optical points is fit¬ 
ted best with a broken power-law—all three extinction curves 
can be used in the modeling of the optical SED. Results are in 
agreement with the X-shooter broadband SED fit and the SMC 
extinction curve. 


The afterglow of GRB 130427A comes with a rich multi¬ 
wavelength data set and has been analysed in detail. It has 
been interpreted either as a combination of a forward and re¬ 
verse shock afterglow contributions rising from a wind ( [Laskar[ 
et al.[[2013] Perley et al.[[2014[ Panaitescm et al.[[2013| l or ISM 


circumburst medium (Maselli et al. [2014 1 . Regardless of the 
interpretation the cooling frequency lies at vc > Vx at the X- 
shooter epoch. Therefore the break we observe in the SED is 
not the cooling break, but may be due to the contribution of both 
forward and reverse shock to the emission ([Perley et al.|2014|l. 


3.8. GRB130603B 


This is the only short GRB in our sample. Photometric data 
used for normalization and photometric SED analysis were ob¬ 
tained from [de Ugarte Postigo et aL] ( [2014| l. The X-shooter SED 
alone can be fitted by a broken power-law model and all three 
extinction curves. Broadband SED is modeled well with a bro¬ 
ken power-law and all three extinction curves. Visual inspec¬ 
tion shows the LMC and MW-type curves do not reproduce the 
blue part of the SED that well, therefore we prefer the SMC- 
type curve. A single power-law fails in reproducing the SED. In 
this case we included the X-band photometric observation in the 
fitting procedure; without this point the fitted pre-break slope 
would be too steep. However, we note that the NIR SED part 
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Fig. 2. SEDs and best-fitted models of the sample. For details, see caption in Figure[^ 
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Fig. 3. SEDs and best-fitted models of the sample. For details, see caption in Figure[^ 


Article number, page 9 of 


17 




















































Table 3. Overview of the quality of fit results 


GRB 

SMC 

EMC 

MW 

Best model 







Phot 

Spec 

100219A 

3.8/10 

41.7/30 

4.3/10 

34.8/30 

6.9/10 

44.2/30 

MW 

EMC 

100418A 

11.2/12 

20.8/23 

10.7/12 

20.2/23 

11.1/12 

20.0/23 

Any 

SMC+ 

100814A 

48.1/33 

70.8/66 

47.8/33 

71.0/66 

47.9/33 

257/67 

Any 

SMC+ 

100901A 

15.4/26 

44.2/41 

14.8/26 

160/47 

14.6/26 

355/47 

Any 

SMC 

120119A 

59.7/47 

194.1/81 

57.5/47 

106.0/81 

79.5/47 

1023/81 

SMC/LMC 

EMC 

120815A 

21.3/21 

26.0/47 

22.2/21 

122.9/47 

20.0/22 

353.1/27 

Any 

SMC 

130427A 

62.6/64 

129.3/147 

62.5/64 

130.0/147 

62.5/64 

123.9/147 

Any 

SMCt 

130603B 

11.5/9 

21.3/23 

11.1/9 

20.9/23 

10.6/9 

21.2/23 

Any 

SMCt 

130606A(“> 

15.9/18 

48.8/31 

16.0/18 

48.8/31 

16.0/18 

48.8/31 

Any 

Any 


Notes. Results of the fits to the broadband SEDs where photometry (phot) and X-shooter spectrum (spec) is used to build the NIR-to-UV part. 
Here we focus on the extinction curve component of the models. For each model we report the ;^^/dof; here we consider only models where 
and Nhi are left as free parameters (see Table [ATT] for details and other models). The last two columns show whether a certain extinction curve 
is found to be preferential over the other two. In the X-shooter case, the better coverage of the blu e par t of the optical SED sometimes allows to 
distinguish between models even if their fits are of a similar quality according to thex^ (see Section [43| : such models are marked with f. 

(a) The case of GRB 130606A is found to have Ay ~ 0. 


cannot be modeled very well (see the residual plot in Figure |^. 
The difference in spectral slopes is A/3 ~ 0.5; fixing the spectral 
difference to = 0.5 therefore does not change the results. A 
similar analysis has been done by |de Ugarte Postigo et al.| ( [2014| l. 
They also find a high extinction of Ay ~ 0.9, albeit a bit lower 
than we do (Ay ~ 1.2). They find a position of the spectral 
break to be near 10*® Hz, while our best fit prefers a value of 
ss 0.8 X 10*® Hz. 


The broadband SED with photometric optical points is fitted 
best with a broken power-law—all three extinction curves can be 
used in the modeling of the optical SED, though the MW curve 
seems to provide a slightly better fit than the other two curves. 
As in the X-shooter broadband fit, the spectral break is found 
in the optical region. However, best-modeled extinction Ay is 
found to be significantly lower. 


Optical afterglow observations of the GRB 130603B are 
sparse. Optical and X-ray light curves appear to evolve achro¬ 
matically after ~ 0.25 days, but are very different prior to that 
time ( |de Ugarte Postigo et al.|20T4| ). During this time a gradual 
steepening from a ^ 1.4 to a ^ 2.4 in X-rays implies an occur¬ 
rence of a jet break ( |Eong et ar]|2014| l. The complicated light 
curve behavior prevents us to use the closure relations. We note 
that both |Eong k al. ( 2014| l and de Ugarte Postigo et al. ( |2014| l 
find a spectral break between the optical and X-ray region, but 
at higher frequencies than we do. 


3.9. GRB130606A 


Photometric data used for normalization and photometric SED 
analysis were obtained from Castro-Tirado et al. P013| l and 
|Afonso et al.| ( |2013] l. The fit to the X-shooter part of the SED 
alone is not very well constrained. Broadband SED is well 
represented by a power-law. No host extinction is necessary 
(Av < 0.01). This is a high-redshift GRB (z = 5.913), there¬ 
fore we had to include the Lya absorption to constrain the UV 
slope. The value of the logAni = 19.9 + 0.3 that we obtain in 
the fit agrees with the one derived by|Hartoog et al.|(|2014|l (i.e., 
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logAffl = 19.94 + 0.01). Eixing the value to logAni = 19.94 
therefore does not change the results. 

Eor the photometric SED we use only filters not affected 
by host’s Lya and intergalactic medium absorption. Broadband 
SED with photometric optical points is fitted best with a single 
power-law and negligible extinction (while upper limits are quite 
high, clearly no extinction is necessary to model the SED). The 
results are in agreement with those obtained from the X-shooter 
broadband SED. 

At the X-shooter epoch the light curve of GRB 130606A 
decays achromatically with a ~ 1.9 both at optical and X-ray 
wavelengths ( Castro-Tirado et al.|2013| l. Our spectral slope and 
temporal decay index are consistent within the case of vc > vx 
and wind circumburst environment. 


4. Discussion 

4.1. Spectroscopic vs photometric SEDs 

To show the advantage of using spectroscmic SEDs in dust anal¬ 
ysis we also modeled photometric SED^ using the photomet¬ 
ric data shown in Eigures[^-[^(but adding X-band magnitudes, 
where available). Due to a poor sampling in blue SED parts, we 
did not use photometric points contaminated by Lya absorption. 
While the difference in the spectral slopes A/3 in the case of pho¬ 
tometric SED studies is usually fixed to A/3 - 0.5, we leave it as 
a free parameter, because the same prescription was used in the 
spectroscopic SED study and therefore the comparison between 
the two analyses is more genuine. The results are reported in 
Table lAA] 

We find that the results of the modeling of photometric SEDs 
are similar in spectral shape to the results of modeling the spec¬ 
troscopic SEDs. However, photometric SEDs can usually be 
modeled by more than one extinction curve; in most cases we 
cannot strongly favor one of the models over the other two. Sim¬ 
ilar conclusions have been reached in other sample studies (e.g., 

^If not stated otherwise, we use the terms “spectroscopic" and “pho¬ 
tometric" to refer to the broadband SEDs which include X-ray data. 
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Fig. 4. Comparison of best-fit parameters of extinction Ay (left) and optical spectral index Po (right) obtained by fitting the broadband SED with 
photometric measurements and X-shooter spectra. Dashed lines represent relations Ay.xshoot = Ay.phot and jSo.xshoot = jSo.phot- Filled and empty 
symbols represent the cases where the broadband SED is best described by a single or broken power-law, respectively. 


Covino et al.|2013 1. The poor resolution of multicolor photom¬ 


etry is not sufficient to detect the smooth features in the SED 
which separate one extinction curve from another. This is in con¬ 
trast to what is found in X-shooter SED analysis. In four cases of 
our sample one extinction curve is strongly preferred according 
to the;^^ statistics. In four additional cases, for which the best fit 
models with different extinction curves are of a similar quality, 
the better coverage of the blue part of the X-shooter SED allows 
us to distinguish the best one by visually inspecting the blue part 
of the SED, which is the most sensitive to the modeling. In these 
four cases the SMC-type extinction curve matched the blue data 
well, while for the other two curves either the 2175 A absorp¬ 
tion was overestimated or the bluest SED part was overpredicted 
by the modeling (see Section]^ for details). The ability to dis¬ 
tinguish between the different extinction curves is summarized 
in Table In the spectroscopic analysis we were thus able to 
single out the best model for all cases but GRB 130606A, whose 
line-of-sight lacks a notable extinction in the first place. 


Differences are also observed in the values of the best-fit pa¬ 
rameters. Eigurej^compares the values of extinction Ay and op¬ 
tical spectral index jSo obtained from both types of analysis. The 
extinction values seem to be systematically larger in the spectro¬ 
scopic analysis, while optical spectral indices do not show any 
preferential deviation. Eixing the difference in spectral slope to 

= 0.5, as it is usually done in these types of analyses (e.g., 
[Greiner et al.|201 f] jCovino et al.|2013l ), does not reduce the dif¬ 
ferences significantly. Due to rather large errors in the values 
derived from the photometric analysis and small sample we can¬ 
not draw strong conclusions about possible trends in deviation. 
However, we emphasize that the photometric SEDs used in the 
analysis cover a very broad wavelength region. The analysis of 
SEDs with less covered spectral range, which are still often used 
in SED modeling, would result in even greater differences from 
the spectroscopic analysis. 

The equivalent hydrogen column densities were not con¬ 
strained very well by fitting only the X-ray part of the SED (see 
Table [A. l) . The reason for this is the low signal-to-noise of the 
X-ray spectrum: to minimize the error due to possible spectral 
evolution and uncertainties in temporal extrapolation, the X-ray 
SEDs were built from rather narrow time intervals. The broad¬ 
band fit resulted in more constrained values of Ah.x- We note 


that the spectral slopes in the X-ray are systematically steeper 
(although still within the error) in the broadband with respect to 
the X-ray-only fits. On the other hand, the fit of the X-shooter 
SED alone usually turned out very bad—only in three cases was 
such a fit both statistically acceptable (i.e., see Section and 
resulted in similar parameter values as in the case of the broad¬ 
band fit. Thus, independent of the quality (resolution) of the 
optical-to-NIR data, the SED modeling is secure only when the 
fit is performed on a broader band SED. 


4.2. Extinction and equivaient hydrogen coiumn densities 


The values of extinction (0 < Ay < 1.2) and equivalent hydrogen 
column density (0.1 < Nu,x [lO^^cm^^j < 2; not including upper 
limits) that we find for our sample are similar to those found in 
GRB lines-of-sight (e.g., Zafar et al.|201 1 [Greiner et al.|201 T 


Watson & Jakobsson|2012f Campana et al.|2012| Covino et al.j 


2013| l. This is illustrated in Figure |5| where we plot both quan¬ 
tities as a function of redshift and compare them to the Swift 
BAT6 complete GRB sample ( [Salvaterra et al.|[2012| jCovinoj 
jet al.|[20T3| ). Most of the events are found to have low extinc¬ 
tion (Ay < 0.3). GRBs 120119A and 130603B are moderately 
extinguished with Ay ~ 1.1, which according to the complete 
Swift BAT6 sample places them into the top 20-25%. The short 
GRB 130603B has the highest extinction in the sample. Dust-to¬ 
gas ratios of our (long GRB) sample are low with respect to the 


ones of the Local Group, which is a well known result (Stratta 
jet al.|[2004) jSchady et al.|201^ [Zafar et ^[201 Ij l. The ratio 


also appears to increase with redshift-since the extinction does 
not evolve much with redshift (see Figure |^, the evolution of 
the ratio is a consequence of the lack of low Nh,x values at high 
redshifts ( [Watson et al.|[20l3] [Covino et al.[[2013| l. The reason 
for the preference of high Nh,x values at high redshifts is not yet 
clear. [Campana et al.[ ( |2()12| l claims that for high-redshift GRBs 
the absorption by intervening systems in the GRB line-of-sight 
contributes much to the measured Ah,x- Alternatively, [Watson 
[et al.[ ( [2013| l claims that the absorption is intrinsic to the GRB 
environment and that the evolution of the (metal) gas column 
density reflects the evolution of cosmic metallicity. As already 
noted by[de Ugarte Postigo et al.[([2014|l, the dust-to-gas ratio of 
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Fig. 5. Extinction and hydrogen equivalent column density values as a function of redshift for our sample. Long GRBs are plotted with empty 
red squares while the short GRB 130603B is given with a filled red square. Empty blue circles represent values obtained in a study of a complete 
sample of long GRBs ( [Covino et al.|2013) . 


the short GRB 130603B is consistent with the Galactic, indicat¬ 
ing that the explosion site for this event differs from a typical 
long GRB site. 

4.3. Extinction curves 

We find that the dust properties preferred for the GRB sight-lines 
are those of the SMC type, with six afterglow SEDs best htted 
by the corresponding averaged extinction curve. Two are best 
described by the EMC curve and one (GRB 130606A) is found 
to have Ay ~ 0. While the preference for the SMC-type of dust 
has been already observed in the sample studies using photo¬ 
metric SEDs (e.g. Kann et al.|2010[ Greiner et al.|201 f] Covino 
|et al. [2013) 1 and the spectroscopic study of Z11, it is still interest¬ 
ing that the well sampled and broadband X-shooter SEDs can be 
modeled with this average extinction curve that well. The lack 
of the 2175 A bump in all but one event (see below) is not sur¬ 
prising; Z11 found that the events with a notable 2175 A feature 
all have rather high extinction values and that the preference for 
the SMC-type dust can be attributed to the observed GRBs being 
biased towards low extinction lines-of-sight. Indeed, most of our 
events have low measured extinction. If the total-to-selective ex¬ 
tinction Ry in the lines-of-sight in GRB hosts were larger than in 
the three assumed extinction curves, the curves would be flatter 
and the derived extinction higher. Such extinction could occur if 
dust were being destroyed by a GRB and grains of smaller size 
are preferentially destroyed ( jWaxman & Drain^|2000| l. How¬ 
ever, while such dust may have been found in some events ( |Per- | 
ley et al.|2008[ Liang & Li|2010 1, as already emphasized in Sec- 


tion |2.2| we find flat extinction curve to be inadequate to describe 
the dust in all lines-of-sight of our GRB sample. 

GRB 120119A is the only event in our sample with a clear 
2175 A absorption bump (see Eigure [^, the feature obs erved 
also in a few other afterglow spectra (Elfasdottir et al.|2009 


[Prochaska et al.|2009||Zafar et al.|2012 |Eynbo et al.|2014 i. In 

our analysis the EMC template provided the best fit to the data 
for this GRB. Still, it is clear from Eigure]^ that the EMC over¬ 
predicts the strength of the 2175 A bump. The failure to find 
a good model is not completely unexpected, since we are using 
merely average extinction curves. Different lines-of-sight in our 


Galaxy or in the two Magellanic Clouds have different extinction 
properties - there are known lines-of-sight in the Galaxy having 
SMC-type dust and vice versa (e.g., see |Elfasd6ttir et al.|2009| for 
review). The analysis thus calls for a more detailed extinction 
model, like the one introduced by [Eitzpatrick & Massa| ( 2007| l, 
where the strength of the bump as well as the UV extinction 
slope are free parameters of the model. The use of more gen¬ 
eral extinction curves in the X-shooter SED analysis will be pre¬ 
sented in a separate work (Zafar et al. in prep.). GRB 120119A is 
a perfect example with excellent data both in optical and X-ray 
frequency range that clearly shows the power of the X-shooter 
data—a photometric SED can hardly differentiate between the 
three models, while the X-shooter data allow us to extract much 
more detailed information about the dust. 


4.4. Neutrai hydrogen coiumn densities 

Eour GRBs in the sample are at a redshift z > 1.7 for which 
the red wings of the Lyo- absorption line enter the X-shooter 
observational window and influence the SED shape. The mod¬ 
eled hydrogen column density values for these bursts are given 
in Table In two cases our values agree with the ones de¬ 
rived from normalized spectra (i.e., GRB 100219A -|Th6ne et al.| 
[20T3| GRB 130606A - [Hartoog et ai:]|20T^ . Eor the other two 
we derive values that are slightly (GRB 120815A - [Krtihler et al. 
|2013| l or significantly (GRB 120119A - Vreeswijk et al. 2014 
in prep) higher from those in the literature. In the case of the 
GRB 120815A, the Lya line is contaminated by vibrationally ex¬ 
cited H 2 lines that form a continuum around this spectral region 
and cause the line to appear stronger ( |Kruhler et al.|20l3] ). The 
contribution of the Hj absorption has already been subtracted 
from the value of log Ahi = 21.95 + 0.15, reported by |Kruhler| 
|et al.| ( [2OT3) . Not taking the molecular hydrogen into account, 
they measure log Am ~ 22.1, closer to our value. 

There are several possible reasons for the further discrep¬ 
ancy. Eirstly, we normalized the spectra of the four GRB af¬ 
terglows and fitted the Lya lines; our best-fit Ahi match very 
well with the values from the literature. Secondly, our spec¬ 
tra are heavily binned. To check the dependency of the results 
on the bin size, we redid the modeling of the X-shooter SED 
part but with smaller bin widths in the region around the line; 
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Fig. 6. Zoomed blue X-shooter part of the broadband SED fit of 
GRB 120119A. Several models are shown, indicating their failure to 
fully reproduce the data. 


the obtained value of the logAHi does not change significantly. 
We also fixed the logAm to the values measured from the nor¬ 
malized spectra (see Table [ATT] and Section]^. The results did 
not change significantly for GRBs 100219A and 130606A. In 
the case of GRB 120815A the fit statistics was worse, but oth¬ 
erwise the parameter values did not change significantly. The 
case of GRB 120119A remains puzzling. The value of log Ahi = 
23.4 + 0.2 found for GRB 120119A seems to be uru'ealistically 
high. While GRB lines-of-sight are characterized by generally 
high hydrogen column densities (|Fynbo et al.|2009 Thone et al. 


|2013| l, the measured values have never surpassed log Ahi = 23.0. 
Indeed, the column density measured from the normalized spec¬ 
trum is found to be 22.5+0.3 (Vreeswijk et al. in prep), but fixing 
this value in the fitting procedure cannot reproduce the blue part 
of the SED (see Figure]^. The main reason for the discrepancy 
is probably due to a low signal-to-noise of the spectral region in 
which the line lies (i.e., the bluest part of the spectrum). The 
uncertainty of our result is further enhanced by the shape of the 
extinction curves that in this case have proved to be inadequate 
to model the extinction. 


4 . 5 . Standard afterglow model 

While this is not the main goal of the paper, here we present 
a brief discussion of our results as analysed with respect to the 
standard afterglow theory. 

A spectral break is necessary to model the SEDs of six GRBs 
in our sample; GRB 100418A, 100814A, 100901A, 100815A, 
130427A and 130603B. An occurrence of a spectral break in 
the SED is not surprising: analysing a sample of ~ 40 after¬ 
glows, Z11 find that ~ 60% of SEDs have a spectral break lying 
somewhere between the optical and X-ray frequencies. Spectral 
slopes for the six GRBs, together with data of the sample pre¬ 
sented by Z11, are plotted in Figure |7] In four of our cases, the 
difference between pre- and post-break slopes of our GRBs is 
in agreement within errors with theoretically predicted value of 
A/3 - 0.5, confirming that both optical and X-ray emission have 
a synchrotron origin and are produced by the same mechanism - 
similar conclusions were reached by Zll. Nevertheless, we find 
two exceptions to the expected behavior. The SEDs of GRBs 



Fig. 7. Comparison of optical and X-ray spectral slopes for the 
cases where spectral break is found to lie in the analysed spectral range. 
Our values (empty red squares) are compared to the sample of |Zafar| 
|et al.| ( (201 l| l (filled black circles). Dashed line represents the relation 
ySx = /6o + 0.5, corresponding to a theoretical prediction in the case the 
observed break is the cooling break. The data point corresponding to 
GRB 130427A, for which the break probably does not correspond to a 
cooling break, is plotted with filled square. 


100418A and 130427A are found to have a break with much 
shallower A/3 ~ 0.3. The afterglow of the GRB 130427A has 
been interpreted either as a combination of a forward and reverse 
afterglow contribution rising from a wind (|Laskar et al.||2013[ 
|Perley et al.||2014[ |Panaitescu et al.||2013| l or ISM circumburst 
medium ( Maselli et al.||2014) l. Regardless of the interpretation, 
the cooling frequency is found to lie above the X-ray frequency 
at the X-shooter epoch, suggesting that the break in the SED is 
not due to the passing of the cooling frequency through the opti¬ 
cal SED. This is supported by the fact that the measured change 
in spectral slopes differs from the one predicted by the standard 
theory for the cooling break. The break we find may thus be due 
to the contribution of both forward and reverse shock to the emis¬ 
sion ( [Perley et al.|2014| l. The case of GRB 100418A is harder to 
explain. Its optical afterglow is characterized by a strong bump 
which is not seen in the X-rays - our spectrum is taken when the 
bump is still strong. A contribution of multiple components in 
the optical may be the cause of the shallow A/3 = 0.31 + 0.09. 
Even greater discrepancy was found for GRB 0802 lOA by Z11 
and De Cia et al. ( |201 l| l, for which the spectral difference is 
A/3 ~ 1. However we recall that we used a sharp spectral break 
in the fitting procedure (see Section [23] l. If the actual spectral 
transition is somewhat milder, occurring over several orders of 
magnitudes, the modeled spectral difference would be underes¬ 
timated by using a sharp break. We also point to rather large 
uncertainties of optical spectral indices. This is because spectral 
breaks lie within or near the optical spectral region and conse¬ 
quently the spectral range needed to constrain the optical slopes 
is quite narrow. 

We briefly checked the available afterglow light curves of 
the GRBs in our sample to see whether the relations between 
spectral and temporal indices are in agreement with theoretical 
expectations in the case of different physical models (e.g.. Table 
1 in Racusin et al.|2009) l. A detailed account on the consistency 
with closure relations for each burst is given in Section]^ While 
some cases can indeed be explained within one of the models’ 
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predictions (e.g., 100814A and 130606A), most of the events 
are not consistent with closure relations. This is a known prob¬ 
lem observed in other sample studies (e.g., Z11, [Zaninoni et al. 
|2013| l and probably reflects our ignorance of the hydrodynamic 
evolution of the outflow and the behavior of the central engine. 


5. Conclusions 


We have presented a detailed SED analysis of a sample of GRB 
afterglows, observed with the VLT/X-shooter spectrograph. Our 
aim was to model broadband afterglow SEDs, composed of the 
NlR-to-ultraviolet X-shooter spectra and SwiftfKKI X-ray ob¬ 
servations, to infer dust properties along the lines-of-sight in 
GRB host galaxies. Only spectra with trustworthy flux calibra¬ 
tion (i.e., for which the spectral slope matches with the one built 
from photometric measurements) were used in the analysis. The 
sample consists of nine GRBs: eight GRBs belong to the long 
and one to the short class. 

The values of extinction (0 < Ay < 1.2) and equivalent hy¬ 
drogen column density (0.1 < A^h.x [lO^^cm^^j < 2) which we 
measure are commonly found in GRB sightlines (e.g., [Cc 


et al.|201? I. Six GRBs in the sample require a spectral break in 
the modeled spectral region. Interestingly, four breaks occur in 
the region which is covered by X-shooter. Most of the events are 
best described by the SMC-type of dust. 

We redid the modeling using the broadband SEDs with the 
NIR/UV photometric measurements in order to see whether 
there are any differences in using the two data sets in the analy¬ 
sis. The derived values of the extinction Ay and spectral slope, 
obtained through spectroscopic and photometric SED analysis, 
for individual events can differ significantly, though no apparent 
trend in the differences is observed. More importantly, the spec¬ 
troscopic data, especially with their excellent coverage of the 
blue SED part, can help us to constrain the dust properties (e.g., 
the extinction curve) much better than the photometric measure¬ 
ments. In addition, we show that independently of the quality 
(resolution) of the optical-to-NIR data, the SED modeling is se¬ 
cure only when the fit is performed on a SED covering broader 
spectral region (in our case, including the X-rays). 

We have demonstrated that spectra, obtained with X-shooter, 
can be successfully used for a detailed SED analysis. X-shooter 
spectra with good flux calibration can be used to constrain the 
extinction curves and can therefore improve our understanding 
of dust in the high redshift Universe. The study will be ex¬ 
panded, when more photometric observations of afterglows be¬ 
come available in the future. 
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